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PURPOSE:  This  technical  note  examines  dissolved  oxygen  (DO)  levels  and  changes  in  river 
microbiology  during  winter,  low  temperatures,  and  periods  of  ice  cover  with  the  objective  of 
providing  guidance  for  winter  water  quality  modeling.  The  results  of  this  study  provide  information 
on  cold-water  aquatic  microbial  activity  and  the  related  coefficients  necessary  to  model  such 
activity.  The  product  of  this  research  can  be  applied  immediately  to  improve  winter  modeling  using 
existing  water  quality  models  such  as  the  Corps’  CE-QUAL-W2. 

BACKGROUND:  Adequate  levels  of  DO  are  critical  to  sustain  aquatic  life  and  maintain  acceptable 
water  quality.  The  most  critical  period  for  dissolved  oxygen  deficits  is  often  assumed  to  be  summer, 
because  low  flows  often  coincide  with  high  water  temperatures.  However,  previous  research  has 
documented  the  occurrence  of  wintertime  oxygen  deficits  in  ice-covered  rivers.  DO  deficits  have 
been  known  to  coincide  with  ice  cover  formation  (Schreier  et  al.  1980;  Whitfield  and  M‘^Naughton 
1986)  and  DO  deficits  can  increase  until  ice  cover  breakup  (Shallock  and  Lotspeich  1974).  Rapid 
decreases  in  winter  DO  have  been  observed  at  river-reservoir  confluences  (James  et  al.  1992, 
McBean  et  al.  1979),  with  large  oxygen  deficits  reported  in  deep  stratified  reservoirs  (Gamble  1971). 

Despite  these  observations,  dissolved  oxygen  processes  in  ice-covered  rivers  have  received  little 
attention.  Given  that  the  solubility  of  oxygen  in  water  is  higher  at  lower  water  temperatures 
experienced  during  winter,  the  same  magnitude  of  oxygen  depletion  due  to  biochemical  oxidation 
has  less  effect  during  the  winter  than  during  the  summer.  In  addition,  the  microbiological  processes 
involved  in  the  consumption  of  oxygen  through  the  decay  of  organic  matter  are  thought  to  be 
insignificant  during  winter  since  they  are  temperature-dependent,  with  slower  rates  at  lower 
temperatures.  Reaeration  during  winter  is  often  assumed  to  be  at  a  minimum  because  ice  covers 
might  provide  a  barrier  to  reaeration.  They  are  also  thought  to  prevent  or  minimize  photosynthesis 
through  the  absorption  or  reflection  of  solar  radiation  in  the  form  of  visible  light. 

However,  low-temperature  microbial  activity  that  occurs  for  extended  periods  in  ice-covered  rivers 
and  reservoirs  can  have  a  significant  impact  on  water  quality  due  to  decreases  in  DO  levels  (Reitner 
et  al.  1997,  James  et  al.  1992,  McBean  et  al.  1979).  This  effect  can  be  magnified  in  rivers  that 
experience  low  flow  during  winter.  Changes  in  a  river’ s  microbial  diversity  and  community  structure 
that  may  occur  during  fall,  winter,  and  spring  temperature  cycles  have  not  been  well-characterized. 
Seasonal  microbial  community  structure  changes  might  influence  the  dominant  metabolic  processes 
occurring  in  a  river.  Degraded  water  quality  manifested  as  reduced  DO  can  result  from  shock 
loadings  that  occur  during  periods  when  the  stability  of  a  river’s  microbial  community  is  low  because 
of  the  lack  of  diversity.  Shock  loadings  include  rainfall  or  snowmelt  runoff  containing  high 
biochemical  oxygen  demand  (BOD)  highway  and  airfield  deicers,  as  well  as  unusually  large 
outflows  from  combined  sewer  systems.  Understanding  the  relationships  among  seasonal  river 
changes  and  river  microbial  community  structure  will  lead  to  improved  predictive  capabilities  for 
winter  water  quality. 
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Quantitative  wintertime  DO  modeling  to  examine  winter  effects  is  difficult  because  conventional 
water  quality  models  (e.g.  U.S.  Environmental  Protection  Agency’s  (EPA’s)  QUAL2E  (EPA  1995) 
and  the  Corps  of  Engineers’  CE-QUAL-W2  (Cole  and  Buchak  1995))  are  extremely  complex.  They 
require  calibration  data  that,  while  sometimes  difficult  to  obtain  in  open-water  conditions,  are 
generally  unavailable  for  ice-covered  rivers.  Water  quality  models  also  provide  uncertain  results 
during  the  winter  because  they  neglect  the  effects  of  reaeration  and  deoxygenation  at  the  low 
temperatures  occurring  during  winter.  The  inclusion  of  winter  dissolved  oxygen  processes  in  simpler 
models  will  allow  greater  flexibility  in  analyzing  and  predicting  riverine  water  quality. 


Riverine  microbial  communities  in  winter  may  be  different  than  those  in  summer,  and  these 
differences  may  have  a  greater  effect  on  water  quality  than  would  be  predicted  by  temperature 
differences  alone.  Colder  water  temperatures  would  be  expected  to  reduce  microbial  activity  in 
winter.  An  approximate  estimate  of  activity  of  a  “pure”  biologically  mediated  process  at  lower 
temperature  could  be  obtained  from  the  Arrhenius  equation: 


K  =  Ae 


-Ea/RT 


(1) 


where  A  is  a  constant,  is  the  activation  energy,  R  is  the  universal  gas  constant,  and  T  is  the 
temperature.  For  many  biochemical  reactions,  this  relation  yields  a  doubling  or  halving  of  rates  as 
temperatures  increase  or  decrease,  respectively,  by  10  °C.  In  practice,  microbial  activity  in  rivers 
and  subsequent  effects  on  water  quality  are  assumed  to  be  insignificant  at  temperatures  below  about 
5  °C. 

Other  factors,  such  as  river  flow  rates,  nutrient  loads,  DO  levels,  ice  cover,  increased  surface  area 
due  to  ice  cover  and  frazil  ice,  and  temperature  cycles  are  influenced  by  winter  conditions.  Less  is 
known  about  interactions  among  these  factors  and  microbial  influences  on  water  quality.  These 
winter  phenomena  can  impact  microbial  communities.  Frazil  ice,  for  example,  may  increase  surface 
area  sufficiently  to  have  a  significant  effect  on  biofilm  development  and  this  may  alter  water  quality 
in  winter  conditions.  Winter  conditions  may  result  in  different  microbial  communities,  which  in 
turn  may  result  in  changes  in  the  dominant  metabolic  processes  in  a  river.  Sucb  changes  may  affect 
winter  river  water  quality  in  addition  to  temperature-driven  rate  changes. 

In  previous  studies  using  soil  samples,  we  have  shown  that  the  effects  of  temperature  cycling  on 
microbial  activity  are  not  intuitive  (Reynolds  et  al.  1998).  Using  a  diurnal  cycle  with  an  amplitude 
of  20  °C,  we  observed  no  significant  difference  in  cumulative  CO2  evolution  between  constant  and 
cyclic  diurnal  treatments  when  soil  temperatures  were  always  above  freezing.  For  the  soil  and  the 
conditions  we  imposed,  using  the  mean  daily  soil  temperature  to  predict  CO2  evolution  appears  to 
be  appropriate  when  temperatures  are  above  freezing. 

However,  when  the  minimum  diurnal  temperature  reached  below  freezing,  even  though  the  mean 
temperature  was  above  freezing  and  the  amplitude  in  the  diurnal  cycling  again  was  20  “C,  CO2 
evolution  increased  in  the  cyclic  diurnal  treatment  relative  to  the  equivalent  constant  temperature. 
In  such  cases,  CO2  evolution  predictions  based  on  daily  mean  temperatures  may  underestimate  CO2 
evolution. 
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Synthesizing  the  above  observations  resulted  in  the  following  hypotheses:  1)  The  effects  of  cyclic 
temperatures  on  biologically  governed  processes  in  streams  and  rivers  may  be  greater  than  are 
currently  recognized,  2)  Omitting  low-temperature  effects  on  microbial  process  rates  may  result  in 
inaccurate  predictions  of  water  quality,  and  3)  Temperature-mediated  changes  may  influence  not 
only  rates,  but  may  also  alter  the  dominant  microbial  community,  thereby  causing  changes  in  the 
primary  metabolic  pathways. 

Unfortunately,  data  that  can  be  used  in  the  development  and  testing  of  hypotheses  related  to  winter 
DO  are  sparse.  The  published  data  that  do  exist  for  a  few  northern  rivers  in  North  America  (including 
the  Athabasca  River  in  Canada  and  the  St.  John  and  Yukon  Rivers  in  the  United  States  and  Canada) 
are  characterized  by  wide  spatial  and  temporal  variations  in  measurement.  Recently,  one-day  winter 
(January  and  March)  DO  and  photosynthetically  active  radiation  (PAR)  data  have  been  reported  by 
Guasch  et  al.  (1998)  for  two  rivers  in  Spain;  however,  no  ice  cover  was  present.  White  and  Melloh 
(1999)  report  high-temporal  resolution  data  for  a  small  river  in  Vermont,  but  the  relatively  pristine 
conditions  there  do  not  yield  data  useful  in  modeling  more  nutrient-rich  river  systems. 

This  study  provided  high  temporal-resolution  water  quality  data  for  an  intermittently  ice-covered, 
relatively  nutrient-rich  river  in  Vermont.  Field  data  and  laboratory  studies  were  performed  to 
evaluate  reaeration  and  deoxygenation  coefficients,  and  microbial  diversity  was  examined. 

METHODS: 

Field:  Hypotheses  were  tested  using  a  combination  of  field  and  laboratory  experiments.  Field  data 
were  collected  from  late  fall  1999  through  spring  2000  on  the  LaPlatte  River  at  Hinesburg,  VT.  The 
LaPlatte  River  watershed  is  53  square  miles,  of  which  47.5  percent  is  agricultural,  40  percent 
forested,  and  8  percent  urban.  Elevations  range  from  1600  ft  mean  sea  level  (MSL)  near  the 
headwaters  in  Hinesburg  to  about  310  ft  MSL  near  the  Hinesburg  Wastewater  Treatment  Plant,  to 
about  100  ft  MSL  at  the  confluence  with  Shelburne  Bay  in  Lake  Champlain.  Water  quality  data 
were  collected  immediately  downstream  from  the  outfall  of  the  Hinesburg  Wastewater  Treatment 
Plant,  which  consists  of  aerated  lagoons  in  series  followed  by  a  chlorine  contact  chamber.  The 
watershed  at  the  sampling  point  is  approximately  18  square  miles,  primarily  forested  with  limited 
agricultural,  residential,  and  commercial  development.  At  this  location,  the  river  is  gently  sloped 
and  meandering,  with  silty  clay  substrate.  The  river  is  surrounded  by  agricultural  land  (com  and 
pasture)  near  the  sampling  site.  A  cheese  plant  is  located  approximately  0.6  mile  upstream  and 
effluent  from  the  cheese  plant  enters  the  Hinesburg  Wastewater  Treatment  Plant  after  treatment  in 
oxidation  ditches  followed  by  clarification. 

A  YSI  600X  probe  was  initially  placed  at  the  center  line  of  the  river.  Total  river  depth  was 
approximately  1  m.  Floats  were  added  to  the  sensor  to  keep  it  from  dropping  to  the  bed  during  low 
flows  and  the  resulting  initial  probe  depth  was  approximately  0.3  m.  The  depth  of  the  probe  varied 
with  flow  during  the  course  of  the  sampling  period.  The  probe  included  sensors  measuring  DO  (mg/1 
and  percent  saturation),  water  temperature  (°C),  conductivity  (mS/cm),  and  pH.  The  probe  was 
calibrated  attbe  Cold  Regions  Research  and  Engineering  Laboratory  (CRREL)  prior  to  installation. 
DO  at  the  time  of  calibration  was  saturated  at  the  local  temperature  and  pressure.  A  three-point  pH 
calibration  was  made  using  buffers  of  4.01,  7.0,  and  10.00.  Specific  conductivity  was  calibrated 
using  a  1 .0-mS/cm  solution.  Post-calibration  was  performed  in  the  field  after  breakup.  A  factory- 
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calibrated  LICOR  spherical  quantum  sensor  (LI-193SA)  measuring  PAR  (iimol/m^  s)  was  placed 
at  the  same  location. 


Both  sensors  were  connected  via 
cable  to  a  Campbell  CRIOX  data 
acquisition  system  located  on  the 
riverbank  (Figure  1 ).  The  data  were 
collected  at  10-min  time  intervals 
for  storage  in  the  data  acquisition 
system.  Stored  data  were  accessed 
and  downloaded  by  computer  via  a 
modem  and  cellular  phone  con¬ 
nected  to  the  data  acquisition  system 
(Figure  2)  at  intervals  ranging  from 
daily  to  every  seven  days.  A  solar 
cell  and  battery  provided  power  for 
the  system.  Air  temperature  and 
battery  voltage  were  also  moni¬ 
tored.  Ice  conditions  (including 
depth  of  snow  on  ice,  ice  morphol- 


Figure  1 .  Data  acquisition  system  and  solar  cell  on  bank  of 
LaPlatte  River,  Hinesburg,  VT 


ogy,  and  bottom  surface  structure) 
were  monitored  periodically.  The  data  were  used 
to  explore  the  effects  of  coefficients  of  reaeration 
on  dissolved  oxygen  during  the  winter  period. 

At  intervals,  water  samples  were  obtained,  stored 
in  sterile  gallon  bottles,  and  transported  to  CRREL 
for  further  testing.  Microbial  samples  were  also 
collected  by  suspending  a  glass-bead  packed  tube 
in  the  water  during  open-water  periods  (Figures  3 
and  4).  The  purpose  of  this  sampling  approach  was 
to  obtain  a  biofilm  on  the  glass  beads  for  charac¬ 
terization,  in  addition  to  the  open-water  samples. 
The  mierobial  samplers  were  autoclaved  prior  to 
sampling.  Microbial  samples  were  obtained  by 
aseptically  rinsing  microorganisms  from  the 
beaded  columns  with  a  known  volume  of  sterile 


water  into  a  sterile  vessel. 


Figure  2.  Data  acquisition  system  showing 


(clockwise  from  top)  Campbell 


Laboratory:  The  laboratory  portion  of  the  study 
used  CRREL’s  temperature-controlled  respirome¬ 
ter  (Figure  5).  Respiration  measurements  using 


CRIOX,  cell  phone,  battery,  and 
modem,  with  voltage  regulator  in  the 
center 


evolved  carbon  dioxide  (CO2)  are  an  estimate  of  general,  non-photosynthetic,  microbial  activity. 


The  respirometer  was  modified  to  allow  for  simultaneous  monitoring  of  samples  incubating  at 


different  temperature  regimes.  Each  temperature  regime  can  be  cycled  in  a  defined  amplitude  and 
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Figure  3.  Flow-through  sampler  designed  to  collect  microbial  samples.  Rubber  stoppers  are  pulled  back 
when  submersed  to  allow  flow  through  and  permit  deposition  on  the  surface  provided  by  the 
glass  beads 
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Figure  4.  Flow-through  microbial  sampler  about  to  be  suspended  in  open  water  upstream  from  the  water 
quality  sensors.  Note  sheet  ice  growth  in  river;  frazil  ice  was  also  present  in  deposits  beneath 
the  sheet  ice 

frequency  schedule  or  held  constant.  Samples  are  incubated  under  low-light  conditions.  The 
respirometer  allows  monitoring  activity  in  near  real  time  by  taking  CO2  headspace  samples 
approximately  every  4  hr.  Results  from  a  range  of  constant  temperature  incubations  may  be  useful 
for  predicting  fall,  winter,  and  spring  processes  for  a  river.  More  commonly,  the  cumulative  CO2 
evolution  curves  and  existing  models  could  be  used  to  describe  microbial  activity.  That  is,  respiration 
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at  a  single  temperature  would  be  used  in  conjunc¬ 
tion  with  Arrhenius  relationships  (e.g.,  Equation  1) 
to  predict  respiration  at  other  temperatures. 

Respiration  tests  were  conducted  to  estimate  if 
river  microbial  community  structure  varied  during 
the  winter.  The  tests  used  water  samples  collected 
for  intermittent  and  permanent  ice  cover  condi¬ 
tions  (6  January  and  15  January  2000)  from  the 
LaPlatte  River.  CO2  evolution  at  ]°C  and  25°C 
was  continuously  monitored.  Appropriate  sterile 
controls  were  included  so  fluctuations  caused  by 
temperature-dependent  CO2  solubility  in  water 
could  be  subtracted. 

In  addition  to  general  heterotrophic  activity  esti¬ 
mated  from  the  CO2  evolution  measurements  that 
are  described  above,  microbial  diversity  was  char¬ 
acterized  and  compared  by  a  modified  fatty  acid 
methyl  ester  (FAME)  procedure  using  whole  cells. 

This  approach  indicates  bacterial  community 
structure  changes  that  could  be  used  in  estimating 
changes  in  the  bacterial  diversity  of  the  river  system. 

The  community  structure  in  the  LaPlatte  River  was 
compared  on  three  sampling  dates  (February  25  and  29  and  April  24)  by  a  fatty-acid-based  technique 
that  uses  culturable  microorganisms  recovered  from  the  glass-bead  sampling  tubes.  These  sample 
times  represented  three  conditions:  when  ice  cover  is  well-established  (25  February),  immediately 
following  ice  breakup  (February  29),  and  after  ice  out  (24  April).  In  brief,  sampling  tubes  (Figure 
3)  were  brought  to  the  laboratory  and  stored  at  4°C.  For  community  structure  comparisons,  the 
rubber  stoppers  were  removed  from  both  ends  of  the  glass  tube  and  the  tube  ends  were  flame-ster¬ 
ilized.  Using  aseptic  techniques,  5  ml  of  sterile,  de-ionized  water  was  slowly  introduced  into  the 
upper  end  of  the  tube.  This  served  as  a  washing  solution  to  remove  the  microorganisms  from  the 
beads.  The  washing  solution  was  uniformly  spread  onto  the  surface  of  the  beads  in  the  end  of  the 
packed  sample  tube  to  limit  channeling  through  the  glass  beads.  Washing  solution  from  the  samplers 
was  collected  into  a  40-ml,  sterile  vial.  After  the  water  was  collected,  the  40-ml  collection  vial  was 
sealed. 

For  plating,  the  40-ml  vial  was  shaken  to  uniformly  mix  the  sample  and  1.0  ml  of  the  sample  was 
removed  for  making  serial  dilutions.  For  each  serial  dilution,  0.1 -ml  aliquots  were  spread-plated, 
in  triplicate,  onto  prepared,  O.IX  tryptic  soy  agar  (TSA)  plates  and  incubated  at  25°C. 

Plates  were  counted  to  estimate  culturable  microorganisms.  For  FAME  characterization,  the 
procedure  used  by  Microbial  Identification  Instrument  (Newark,  DE)  (MIDI)  (Sasser  and  Wichman 
1991).  The  standard  procedure  is  a  five-step  process: 


Figure  5.  CRREL’s  temperature-controlled 
respirometer 
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•  Cell  harvesting:  Cells  are  harvested  to  yield  sufficient  fatty  acids  to  be  readily  measurable. 

•  Saponification;  Cells  are  subjected  to  a  strong  methanolic  base  and  heat.  This  kills  and 
lyses  the  cells,  cleaving  fatty  acids  from  the  cell.  The  free  fatty  acids  are  converted  to  their 
sodium  salts. 

•  Methylation:  Fatty  acid  salts  are  esterified  by  adding  a  methyl  group  to  increase  volatility 
for  chromatographic  analysis. 

•  Extraction  and  base  wash:  Free  fatty  acids  and  cellular  debris  are  removed  from  the 
sample. 

•  Analysis:  The  fatty  acid  sample  is  analyzed  by  gas  chromatography  and  the  resulting  fatty 
acid  profile  is  compared  to  a  library. 

For  the  purposes  of  this  study,  the  entire  lawn  of  microorganisms  was  treated  as  a  sample,  rather 
than  individual  isolates.  Principal  component  analysis  (PCA)  techniques  were  then  used  to  compare 
samples. 

At  the  initial  sample  collection  event,  the  effects  of  the  duration  that  the  packed-bead  samplers 
remained  in  the  river  was  evaluated.  The  purpose  of  this  was  to  estimate  the  time  needed  to  develop 
a  biofilm  on  the  beads.  A  control  (no  river  contact),  and  packed-bead  samplers  that  had  been 
submerged  for  0.5  and  3.0  h  were  evaluated. 

PCA  was  used  to  compare  fatty  acid  profiles  among  sample  times.  PCA  is  a  statistical  technique 
that  is  used  to  extract  the  main  relations  in  data  of  high  dimensionality,  as  is  found  in  comparing 
fatty  acid  profiles  obtained  by  gas  chromatography  (GC)  analysis. 

RESULTS  AND  DISCUSSION:  Water  quality  data  were  collected  continuously  at  the  Hinesburg 
site  between  22  December  1999  and  31  March  2000.  Intermittent  ice  covers  occurred  between  24 
and  26  December,  and  on  28  December,  29  December,  31  December,  and  6  January.  A  permanent 
seasonal  ice  cover  began  on  about 
12  January  2000.  By  15  January, 
the  ice  cover  was  10  cm  in  thick¬ 
ness  and  substantial  amounts  of 
frazil  ice  were  deposited  beneath 
the  solid  ice.  Algae  were  visible 
both  within  the  frazil  ice  and  in 
accumulations  at  the  surface  in 
open-water  areas  (Figure  6).  This 
ice  cover  remained  in  place  until 
between  27  and  28  February  2000 
when  the  ice  cover  broke  up  fol¬ 
lowing  a  snowmelt  and  rainfall 
event.  Open  water  was  observed  on 
29  February. 

Water  quality  data.  The  water  qual¬ 
ity  data  for  ice-covered  and  non¬ 
ice-covered  periods  are  presented 


Figure  6.  Algae  visible  at  the  surface  at  the  ice/open  water 

interface,  and  beneath  the  ice  cover  downstream  (see 
arrows  above).  Frazil  ice  was  packed  to  the  bed  along 
the  banks 
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in  Figure  7.  Statistical  analyses  were  performed  using  SAS®  statistical  analysis  software  (SAS® 
Institute  1 996).  The  sample  populations  for  all  of  the  measured  variables  are  statistically  signifi¬ 
cantly  different  (p  <  0.05)  between  the  ice-covered  and  open-water  periods  as  determined  using  the 
Wilcoxon  rank  sums  test.  Water  temperature,  DO,  pH,  and  conductivity  conditions  observed  during 
the  ice-covered  period  are  relatively  stable  compared  to  open-water  conditions  and  are  similar  to 
those  reported  by  White  and  Melloh  (1999).  The  data  demonstrate  a  stabilizing  effect  of  ice  cover 
on  water  quality  parameters  that  may  provide  some  biological  benefits  to  lotic  communities.  This 
idea  was  suggested  by  Hynes  (1970);  i.e.,  ice  covers  provide  some  biological  benefits  to  lotic 
communities  through  stabilization.  Hynes  also  noted  that  scour  induced  during  ice  cover  breakup 


Figure  7.  Data  for  open-water  versus  ice-covered  periods,  22  December  1 999  through  31  March  2000. 

Mean  value  is  indicated  by  +  in  box,  which  depicts  25th  percentile  (lower  line),  75th  percentile 
(upper  line),  and  median  value  (center  line) 
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may  hold  some  hazards  to  riverine  communities;  this  study  shows  that  high  variability  in  tempera¬ 
ture,  DO,  pH,  and  conductivity  is  also  associated  with  breakup. 

Descriptive  statistics  for  the  water  quality  variables  are  provided  in  Table  1 .  Correlation  tests  were 
also  performed  and  during  the  open-water  periods,  water  temperature  was  moderately  correlated 
with  DO  percent  saturation  ir^=  0.53),  pH  (r^=  0.69),  and  air  temperature  (^2=  0.66).  Percent 
saturation  DO  was  moderately  correlated  with  pH  (r^=  0.77).  Conductivity  and  pH  were  also 
correlated  with  depth  (r^=  -0.57  and  -0.70,  respectively). 


Table  1 

Descriptive  Statistics  for  Water  Quaiity  Data 

Open-Water  Period 

Variable 

N 

Mean 

Std  Dev 

Minimum 

Maximum 

Water  Temperature  (X) 

6747 

1.93 

1.83 

-0.02 

8.03 

Conductivity  (mS) 

6747 

0.15 

0.03 

0.06 

0.26 

DO  (Percent  saturation) 

6747 

113.2 

5.96 

89.60 

127.90 

DO  (mg/L) 

6747 

15.68 

0.78 

11.16 

17.57 

Depth  (m) 

6747 

0.23 

0.15 

-0.04 

0.93 

pH 

6747 

8.10 

0.09 

7.84 

8.38 

PAR  (pmol/m^) 

6747 

87.91 

168.05 

0.60 

1731.00 

Air  Temperature  (°C) 

6747 

1.67 

6.23 

-18.54 

19.37 

Ice-Covered  Period 

Water  Temperature  (X) 

7733 

0.01 

0.03 

-0.06 

0.20 

Conductivity  (mS) 

7734 

0.21 

0.04 

0.08 

0.30 

DO  (Percent  saturation) 

7734 

107.54 

2.73 

100.60 

114.60 

DO  (mg/L) 

7734 

15.70 

0.40 

14.66 

16.76 

Depth  (m) 

7734 

0.39 

0.19 

0.21 

1.39 

pH 

7734 

7.97 

0.06 

7.81 

8.10 

PAR  (pmol/m^) 

7734 

10.36 

21.97 

-2.39 

348.10 

Air  Temperature  (‘'C) 

7734 

-8.68 

8.30 

-32.54 

9.66 

The  PAR  measurements  during  open-water  periods  indicate  that  adequate  incident  sunlight  is 
reaching  the  sensor  level  (approximately  1  ft  deep)  to  allow  for  photosynthetic  activity.  Open-water 
mean  integrated  irradiance  was  7.6  mol/m^  day,  with  a  daytime  average  integrated  irradiance  of  16 
mol/m^  day.  The  maximum  observed  PAR  was  1731  [imol/m^  s  on  18  March.  These  levels  are 
comparable  to  open-water  wintertime  values  reported  by  Guasch  et  al.  (1998)  as  well  as  to 
summertime  measurements  for  shaded  rivers. 

During  the  ice-covered  period,  the  mean  integrated  irradiance  was  0.9  mol/m^  day,  with  a  maximum 
observed  PAR  of  348  pmol/m^  s  on  24  December.  Following  the  formation  of  the  permanent 
seasonal  ice  cover  on  12  January,  the  mean  integrated  irradiance  was  0.7  mol/m^  day.  The  daytime 
average  integrated  irradiance  was  1.4  mol/m^  day,  with  a  mean  observed  PAR  of  16.2  jimol/m^  s. 
The  maximum  observed  PAR  was  269  pmol/m^  s  on  13  January,  just  after  ice  cover  formation. 
These  values  are  somewhat  lower  than  those  reported  by  White  and  Melloh  (1999):  mean  daily 
integrated  irradiance  of  3.3  mol/m^  day  and  mean  and  maximum  PAR  of  488  pmol/m^  s,  perhaps 
because  of  the  presence  of  frazil  ice  deposited  beneath  the  ice  cover. 
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During  the  ice-covered  period,  poor  correlations  were  found  among  the  variables  water  temperature, 
PAR,  DO,  pH,  and  conductivity,  except  in  three  cases:  depth  was  found  to  be  moderately  correlated 
with  DO  concentration,  (r^=  0.50),  DO  percent  saturation  (r^=  0.50),  and  pH  (r^=  -0.56). 

In  comparison,  Vincent  et  al.  (1998)  reported  PAR  values  of  about  80  [imol/m^s  at  10  m  (32.8  ft) 
depth  in  Lake  Vanda  and  about  8  iimol/m^  s  at  5  m  (16.4  ft)  depth  in  Lake  Bonney,  both  freshwater 
Antarctic  Lakes  exhibiting  primary  productivity.  Thus,  despite  the  relatively  low  PAR  levels 
observed  in  the  LaPlatte  River  during  the  ice-covered  period,  photosynthetic  activity  could  be 
expected  to  continue.  In  fact,  Marchand  (1996)  indicates  that  some  aquatic  plants  appear  to  adapt 
to  photosynthesis  under  low  light  levels,  and  Guasch  et  al.  (1998)  suggest  the  same  for  other  biota. 

The  range  of  water  temperature  variation  was  small  during  the  ice-covered  period,  less  than  0.2  °C 
per  day.  Supercooling  (water  temperatures  slightly  less  than  0  °C)  was  noted  during  the  December 
and  January  open-water  periods,  consistent  with  observations  of  frazil  ice  beneath  the  ice  cover  near 
the  sampling  site.  Supercooling  was  also  noted  during  the  early  part  of  the  ice-covered  periods, 
indicating  that  turbulent  open-water  areas  existed  upstream  at  the  time. 

The  observed  DO  levels  were  supersaturated  during  the  ice-covered  period,  averaging  15.7  mg/L 
and  107.5  percent  of  saturation.  These  values  compare  favorably  to  the  data  reported  by  White  and 
Melloh  (1999),  which  averaged  16.4  mg/L  and  113.7  percent  of  saturation  during  the  late  winter 
ice-covered  period.  The  generally  high  DO  levels  may  reflect  highly  efficient  reaeration  at  low 
temperatures,  photosynthetic  input,  or  both.  Diurnal  cycling  of  DO  was  observed  throughout  the 
ice-covered  period. 

The  low-amplitude  DO  fluctuations  observed  are  comparable  to  those  reported  by  Simonsen  and 
Harremoes  (1978)  and  White  and  Melloh  (1999).  Fluctuations  in  water  temperature  do  not  fully 
explain  the  observed  cycling,  indicating  a  biological  component.  Typical  diurnal  DO  cycles  during 
the  ice-covered  period  resemble  that  shown  in  Figure  8,  for  20  January  2000.  DO  levels  are  lowest 
during  the  nighttime  hours,  begin  to  rise  through  the  morning  hours,  and  peak  sometime  in  early  to 
mid-afternoon.  As  the  ice-covered  period  progresses,  the  DO  cycle  becomes  less  well-defined 
(Figure  9),  but  the  trend  remains  the  same.  The  December  open-water  curves  resemble  the  typical 
curve  shown  in  Figure  8,  as  do  those  for  April  and  May. 

However,  DO  cycling  during  the  period  following  breakup  (after  28  February)  does  not  follow  this 
trend.  Instead,  DO  reaches  its  lowest  point  at  around  sunset  and  rises  to  peak  in  the  morning 
(Figure  10).  DO  concentration  and  water  temperature  are  also  much  better  correlated  during  this 
period  (r^=  -0.82)  than  during  the  ice-covered  period  (t^=  -0.08),  since  the  barrier  to  heat  transfer 
provided  by  the  ice  is  absent.  Following  ice  cover  breakup,  DO  peaks  earlier  in  the  day  than  during 
the  ice-covered  period.  This  is  an  indication  that  ice  cover  breakup  induces  changes  in  the  system 
that  drastically  reduce  microbial  impacts.  Such  changes  could  include  the  loss  of  microbial 
communities  living  within  the  frazil  ice  deposit  and  the  scouring  or  erosion  of  microbial  commu¬ 
nities  living  on  bed  or  bank  surfaces.  The  DO  cycling  observed  following  breakup  is  quite  similar 
to  that  observed  by  White  and  Melloh  (1999)  in  the  relatively  pristine  Sleeper’s  River,  VT,  where 
biological  activity  was  thought  to  be  minimal. 
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Figure  9.  Less  well-defined  DO  cycling  occurs  during  later  ice-covered  periods,  as  illustrated  for  the  period 
2-4  February  2000 
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Figure  10.  Typical  DO,  PAR,  and  temperature  data  following  ice  cover  breakup  (21-24  March  2000) 


Respirometer  Data.  Cumulative  evolved  CO2  respiration  data  obtained  from  the  respirometer  study 
are  shown  in  Figure  1 1 .  Water  samples  collected  on  6  and  15  January  were  evaluated  at  1  °C  and 
25  °C.  On  6  January,  the  river  had  only  intermittent  ice  cover,  but  by  15  January  the  ice  cover  was 
approximately  10  cm  thick  at  the  sampling  site.  Both  samples  had  similar  respiration  patterns  at  1 
°C.  This  was  shown  by  a  lag  period  of  approximately  100  hr,  followed  by  an  increase  in  respiration 
from  100  hr  to  270  hr  to  a  higher  but  relatively  constant  rate.  At  1  °C,  samples  from  both  6  and  15 
January  behaved  similarly  with  respect  to  CO2  evolution,  each  evolving  a  cumulative  total  of 
approximately  50  p-l  C02/ml  during  the  270-hr  incubation,  and  each  having  a  lag  time  of  approxi¬ 
mately  100  hr. 

Additional  aliquots  of  the  water  samples  collected  on  6  and  15  January  were  also  incubated  at  25 
°C  (Figure  IIF  Using  constant  temperature  baths  on  the  respirometer,  these  incubations  were 
conducted  concurrently  with  the  1  °C  incubations.  Data  for  both  sample  collection  times  again  show 
a  lag  period  of  approximately  1 00  hr.  However,  from  1 00  to  270  hr,  there  was  greater  CO2  evolution 
from  the  water  sample  collected  on  15  January  under  ice  cover  than  from  the  sample  collected  on 
6  January  in  open  water.  Cumulative  CO2  evolution  values  were  approximately  330|ll  C02/ml  and 
170  pi  C02/ml  for  the  15  January  and  6  January  samples,  respectively.  Because  the  CO2  evolution 
response  to  warmer  incubation  temperatures  was  significantly  different  between  samples  collected 
on  6  and  1 5  January,  these  data  strongly  suggested  that  there  were  different  microbial  communities 
in  the  river  during  open  water  and  ice-covered  periods. 
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Figure  1 1 .  Cumulative  respiration  data  for  samples  obtained  6  and  15  January  2000 


A  possible  explanation  for  greater  CO2  evolution  from  the  15  January  early  ice-covered  sample  than 
from  the  6  January  open-water  sample  is  that  a  frazil  ice  matrix  had  formed  by  1 5  January,  providing 
structure  for  biofilm  attachment  and  growth.  Visual  observations  on  15  January  revealed  algae  both 
within  the  frazil  ice  and  in  accumulations  at  the  surface  in  open-water  areas  (Figure  6).  Either  algae 
or  bacteria  present  in  such  a  biofilm  would  supply  greater  carbon  to  the  system.  Because  the 
respirometer  has  insufficient  light  to  support  algal  photosynthetic  processes  in  the  incubation 
vessels,  catabolic  bacterial  processes  would  dominate  during  respiration  studies.  Greater  release  of 
CO2  from  the  ice-cover  sample  taken  on  15  January  would  support  the  idea  of  greater  biomass 
carbon,  whether  algal,  bacterial,  or  both,  in  the  ice-covered  sample  relative  to  the  open-water  sample. 

The  high  DO  measurements  are  relatively  consistent  during  late  ice-covered  periods  (Figure  9) 
compared  to  open  water  or  following  ice  cover  breakup  (Figure  10).  The  DO  data  suggest  either 
greater  DO  sources  or  lesser  DO  sinks  under  an  ice  cover  relative  to  open  water.  Greater  CO2 
evolution  from  the  respirometer  data  (Figure  11)  suggest  greater  carbon  in  the  water  sample  taken 
during  ice  cover,  again  probably  due  to  algal  biomass  contained  within  the  frazil  ice  matrix. 

Microbial  Data.  In  addition  to  the  respiration  tests,  microbial  communities  were  characterized  with 
a  modified  fatty  acid-based  procedure,  FAME.  This  comparison  is  based  on  eulturable  microorgan¬ 
isms  that  were  recovered  from  the  glass-bead  samplers  and  integrates  the  fatty  acid  profiles  of  all 
eulturable  microorganisms  from  a  sample  into  a  single  profile.  In  effect,  this  approach  treats  an 
entire  sample  as  an  “organism”  and  then  compares  the  similarity  or  dissimilarity  among  samples. 
Because  this  approach  requires  growth  on  0.1  X  TSA,  the  initial  data  favors  bacterial  changes  over 
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fungal  or  algal  changes.  These  samples  represented  periods  with  a  well-established  ice  cover  (25 
February),  immediately  following  ice  breakup  (29  February),  and  over  six  weeks  after  ice  out  (24 
April). 

The  PCA  results  for  the  fatty  acid  data  obtained  from  culturable  microorganisms  recovered  using 
the  glass-bead  samplers  are  shown  in  Figure  12.  Ellipses  have  been  added  simply  to  better  show 
groupings.  For  the  purposes  of  this  study,  the  importance  of  these  data  is  that  the  culturable  microbial 
communities  from  the  river  clearly  separate  into  three  distinct  groups  among  the  three  sampling 
times.  These  data  corroborate  the  respirometer  data  and  support  the  concept  that  there  are  seasonal, 
ice-cover-related  microbial  community  structure  differences  in  the  river. 


Figure  12.  Results  of  the  principal  component  analysis  of  fatty  acid  data  obtained  from  culturable 
microorganisms  collected  by  the  glass  bead  samplers  during  the  ice-covered  period, 
immediately  after  ice  cover  breakup,  and  well  after  ice  cover  breakup 

Reaeration  Coefficients.  The  concentration  of  DO  in  a  river  at  any  time  is  a  function  of  the 
deoxygenation  resulting  from  the  biochemical  oxidation  of  organic  matter  in  the  river  and  its 
sediment  and  the  reaeration  provided  by  tributaries,  surface  runoff,  groundwater  inflow,  photosyn¬ 
thesis,  and  atmospheric  absorption  of  oxygen  (see,  e.g.,  Tchobanoglous  and  Schroeder  (1987)). 
Freeze  and  Cherry  (1979)  note  that  few  measurements  of  DO  in  groundwater  exist,  but  that 
groundwater  found  in  sandy  or  gravelly  soils,  or  where  little  soil  overlies  fractured  rock  can  be 
oxygenated,  while  shallow  groundwater  in  silty  or  clayey  soils,  similar  to  those  in  the  project  area, 
is  commonly  deoxygenated.  The  component  of  reaeration  associated  with  atmospheric  absorption 
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depends  on  the  DO  saturation  level,  the  DO  concentration  of  the  water,  the  air-water  DO  concen¬ 
tration  gradient,  the  area  of  the  water  surface,  and  hydraulic  factors  such  as  turbulence  and 
dispersion.  Its  saturation  level  controls  the  maximum  level  of  oxygen  at  the  surface  of  the  water, 
which  is  a  function  of  temperature  and  barometric  pressure.  However,  within  the  water  column,  the 
water  can  become  supersaturated  by  as  much  as  200  percent  (Thomann  and  Mueller  1987)  due  to 
photosynthetic  oxygen  production. 

The  DO  in  a  river  system  is  utilized  in  aerobic  respiration  as  the  terminal  electron  acceptor,  in  the 
oxidation  of  organic  and  nitrogenous  material,  and  in  the  oxygen  demand  exerted  in  the  benthic 
community.  The  opposing  processes  of  photosynthesis  and  respiration  often  result  in  a  typical 
diurnal  cycle  such  as  that  shown  in  Figure  10.  The  level  of  DO  in  the  sediments  is  related  to  redox 
potential  (E,,)  in  sediments  (Fletcher  1979),  which  is  generally  positive  in  the  well-oxygenated 
surface  made  up  of  larger  particles  with  little  organic  matter  (i.e.,  favoring  oxidation  reactions)  and 
negative  in  the  deeper  sediments  (i.e.,  favoring  reduction  reactions)  where  the  particle  size  is  finer 
and  there  is  more  organic  matter,  thereby  affecting  the  chemical  environment  (Figure  5).  Marchand 
(1996)  notes  that  if  DO  drops  to  very  low  levels  (as  in  lake  bottoms  during  winter),  some  organisms 
switch  to  anaerobic  respiration  using  glycolysis. 

Successful  modeling  of  DO  in  rivers  under  open-water  and  warm-weather  conditions  has  been 
performed  using  both  simple  and  complex  models.  However,  it  appears  that  existing  deterministic 
and  stochastic  models  for  ice-covered  rivers  underestimate  reaeration  during  ice-covered  conditions 
(e.g.,  Landine  (1970),  Chambers  et  al.  (1993)),  do  not  adequately  address  the  hydraulic  effects  of 
the  ice  cover,  and  tend  to  neglect  or  minimize  the  importance  of  ecological  effects  such  as  microbial 
succession  during  winter.  The  modeling  of  DO  in  rivers  was  first  proposed  by  Streeter  and  Phelps 
in  their  classic  1925  paper,  which  presents  a  one-dimensional,  steady-state  analytical  model  to 
describe  the  rate  of  change  in  DO  with  time  as  a  function  of  both  biochemical  oxidation  and 
reaeration.  Based  on  field  observations,  Streeter  and  Phelps  assumed  that  the  concentration  of  DO 
in  a  river  at  any  time  is  a  function  of  deoxygenation  and  the  reaeration  provided  by  the  atmospheric 
absorption  of  oxygen.  Changes  in  DO  due  to  photosynthesis,  sediment  oxygen  demand,  and 
respiration  are  neglected,  and  the  biochemical  oxygen  demand  is  assumed  to  decrease  monotonically 
over  time; 


^ 

~~dt~  dt  ~ 


(2) 


where  t  is  time  in  days,  is  the  deoxygenation  in  mg/1,  is  the  coefficient  defining  the  rate  of 
deoxygenation,  and  L  is  the  biochemical  oxygen  demand  (also  known  as  the  ultimate  BOD,  or 
BOD^).  A  first-order  process  is  also  used  to  model  reaeration: 


dD^ 

dt 


=  -K,D 


(3) 
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where  is  the  deoxygenation  in  mg/1  and  is  the  coefficient  defining  the  rate  of  reaeration.  The 
oxygen  saturation  deficit  D  (mg/I)  is  the  sum  of  the  deoxygenation  (D^)  and  the  reaeration  (D)  at 
any  time  t,  so  that  the  rate  of  change  of  the  oxygen  deficit  is  a  linear,  first-order  differential  equation: 


^  =  ^  +  ^  =  KjL-K,D 
dt  dt  dt 


Typically,  Equations  2  and  4  are  rewritten  as  a  set  of  coupled  first-order  equations  that  describe  that 
change  in  oxygen  deficit  with  distance: 


u  —  +  K^L  =  0  (5) 

dx  ^ 


and 


u  —  +  KaL-K,{Cs-C)^^  (6) 

dx 

where  u  is  velocity,  C  is  the  DO  concentration  at  any  point,  and  is  the  DO  concentration  at 
saturation,  which  varies  with  temperature  and  pressure.  The  oxygen  deficit  D  is  equal  to  C^-C. 
Variations  of  these  equations  are  commonly  used  to  calculate  DO  in  rivers  downstream  from 
wastewater  point  sources  (see,  e.g.,  Tchobanoglous  and  Schroeder  (1987),  and  Thomann  and 
Mueller  (1987)).  White  (1998)  utilized  a  simple  finite  element  model  to  represent  Equations  5  and 
6.  This  model  is  modified  in  the  present  study  to  examine  diurnal  cycling  of  dissolved  oxygen  by 
substituting  dx/dt  for  u  so  that  changes  in  the  variables  are  computed  over  time  rather  than  distance. 
Reaeration  coefficients  can  then  be  examined  during  different  time  periods  for  best  fit. 

Most  models  include  additional  sink  and  source  terms  to  account  for  the  impacts  of  photosynthesis, 
respiration,  and  other  physical  processes  affecting  oxygen  levels.  However,  in  practice  these  sink 
and  source  terms  are  often  neglected,  with  the  result  that  all  of  the  physical  processes  affecting 
dissolved  oxygen  are  included  in  the  reaeration  and  deoxygenation  coefficients.  Little  research  has 
been  done  regarding  appropriate  K^and  Kj  values  at  low  temperatures  in  natural  rivers,  and  as  Eheart 
and  Park  ( 1 989)  noted,  the  temperature  adjustment  factors  commonly  used  are  not  well-established. 
Streeter  and  Phelps  (1925)  estimated  and  A"^for  December  through  February  on  the  Ohio  River 
for  water  years  1915  through  1917.  Their  mean  values  for  and  of  0.3  d'*  and  0.05  d"’, 
respectively,  for  two  reaches  during  the  winter  months  were  statistically  significantly  lower  than 
the  summer  values  (White  1998).  McBean  et  al.  (1979)  neglected  ^^when  modeling  the  St.  John 
River  (Maine  and  Canada)  in  February,  and  used  a  value  of  0.26  d"'  in  November.  They  used  a 
value  of  0.064  d'^  during  both  November  and  February.  Pietroniro,  Chambers,  and  Ferguson  (1998) 
used  /ir^.=0.001  d“'  for  ice-covered  reaches  of  the  Athabasca  River,  following  the  work  of  McDonald, 
Holley,  and  Goudey  (1989),  who  neglected  KJor  the  same  river. 
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The  present  study  utilized  the  observed  data  for  the  LaPIatte  River  to  determine  coefficients  of 
reaeration  and  deoxygenation  during  the  ice-covered  period.  The  diurnal  cycling  of  the  data  clearly 
showed  some  relationship  between  DO  and  PAR,  and  hence  it  seemed  reasonable  to  include  the 
effects  of  PAR  when  computing  DO.  Rather  than  add  a  source  or  sink  term,  it  seemed  practical  to 
modify  the  A'^term,  since  changes  in  the  level  of  reaeration  would  be  the  primary  result  of  changes 
in  PAR.  A  simple  linear  relationship  was  used  to  account  for  these  effects: 

=  (7) 


where  K  ^  is  the  reaeration  coefficient  at  time  t  modified  to  reflect  PAR  and  a  is  a  coefficient 
describinrthe  relationship  between  PAR  at  time  t-l  (P,, )  where  I  is  a  specified  lag  time  ranging 
from  0  for  open-water  conditions  to  about  3  hr  for  late  ice-covered  period  conditions.  Because 
continuous  BOD  data  were  not  available,  BOD  was  assumed  to  be  constant.  Initial  values  of  = 
0.3  d'^  and  =  0.06  d'^  were  selected.  By  trial  and  error,  the  smallest  least-squared  error  during 
ice-covered  periods  was  achieved  using  =  0.7  d'^  and  =  0.1  d'^,  with  cl  =  0.02  and  a  lag  time 
of  0.12  day.  Lag  time  and  decreased  for  wintertime  open-water  periods,  while  increased 
slightly. 

Despite  the  simplicity  of  the  model  and  its  assumptions,  computed  DO  is  quite  close  to  observed 
DO,  and  does  exhibit  diurnal  cycling  similar  to  that  observed.  For  example,  during  the  open-water 
period  of  1  to  6  January  2000,  DO  varied  between  about  1 3. 1  mg/L  and  14.2  mg/L  with  clear  diurnal 
variations,  and  averaged  13.7  mg/L.  For  comparison,  the  computed  DO  during  the  same  period 
ranged  from  13.4  to  14.2  mg/L  with  a  mean  value  of  13.7  mg/L.  The  computed  and  observed  DO 
values  for  this  period  are  shown  in  Figure  13.  Here,  =  0.5  d"'  and  Kj=^0A2  d'^  with  a  =  0.02 
and  no  lag  time.  If  the  relationship  shown  in  Equation  7  is  not  used,  becomes  constant  and  diurnal 
cycling  is  not  modeled  (Figure  13). 

During  the  ice-covered  period,  smaller  diurnal  fluctuations  in  DO  were  observed.  Typical  model 
output  for  the  later  ice-covered  period,  31  January  to  15  February  2000,  is  shown  in  Figure  14,  with 
=  0.7  d'^  Kj  =  0A  d'*,  a  =  0.02,  and  I  =  0.12  day.  During  this  period,  DO  ranged  from  13.0  mg/L 
to"^  14.0  mg/L,  with  a  mean  value  of  1 3.6  mg/L.  The  low  levels  of  observed  PAR  prior  to  9  February 
(daily  maximums  less  than  15  iimol/m^  s)  induce  small  yet  discernible  changes  in  computed  diurnal 
cycling  of  DO.  Increased  levels  of  PAR  on  10, 12, 13,  and  15  February  result  in  increased  levels  of 
computed  DO  for  those  days.  The  relatively  large  increases  in  observed  DO  on  8  and  13-14  February 
are  more  closely  correlated  with  decreases  in  air  temperature  than  changes  in  PAR,  and  thus  the 
model  slightly  underestimates  DO  on  these  days.  Computed  DO  for  the  period  31  January  to  15 
February  2000  ranged  from  13.1  mg/L  to  13.8  mg/L  with  a  mean  value  of  13.6  mg/L.  However,  the 
error  is  much  less  than  if  constant  =  0.7  d"’  is  used  (Figure  14),  which  ranged  from  13.1  mg/L 
to  14.3  mg/L  with  a  mean  value  of  13.9  mg/L. 

CONCLUSIONS  AND  RECOMMENDATIONS:  The  purpose  of  this  study  was  to  investigate 
riverine  dissolved  oxygen  concentration  during  winter,  including  the  ice-covered  period,  with 
emphasis  on  microbial  processes  contributing  to  reaeration  and  deoxygenation.  Data  showed  that 
DO  patterns  in  a  small  Vermont  river  are  different  during  an  ice-covered  period  than  during 
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Figure  13.  Observed  versus  computed  DO  for  the  open-water  period  1-6  January  2000.  For  comparisi 
DO  computed  using  constant  Kr=  0.5  cf’  is  also  shown.  Observed  PAR  is  shown  for  the  s( 
period 
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Figure  14.  Observed  versus  computed  DO  and  observed  PAR  and  air  temperature  during  the  late 
Ice-covered  period  1-15  February  2000.  For  comparison,  DO  computed  using  constant 
Kr=  0.7  is  also  shown 
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open-water  periods.  In  addition,  diurnal  cycling  of  DO  persists  during  the  ice-covered  period,  and 
appears  to  be  related  to  PAR  in  a  manner  that  suggests  photosynthetic  activity  may  be  contributing 
to  reaeration.  Visual  observations  at  the  time  of  ice  cover  formation  showed  algae  present  in  the  ice 
cover  and  within  the  frazil  matrix  trapped  beneath  the  ice  cover.  It  is  possible  that  the  porous  frazil 
structure  provides  a  surface  for  biofilm  formation  that  supports  microbial  communities  in  a  different 
manner  than  the  typical  lotic  environment.  Microbial  communities  were  found  by  two  independent 
measures  to  be  different  under  an  ice  cover  than  in  open  water. 

Previous  modelers  have  either  neglected  reaeration  during  winter  or  have  used  very  small  values 
for  the  reaeration  coefficient  K^.  Data  showed  DO  levels  to  be  at  or  above  saturation  levels  for  the 
ice-covered  period,  indicating  that  reaeration  could  be  important.  A  simple  finite  element  model 
based  on  the  classic  Streeter-Phelps  model  was  applied  to  the  field  data  to  evaluate  reaeration  and 
deoxygenation  coefficients.  The  use  of  a  simple  linear  relationship  in  the  model  between  PAR  and 
K  to  account  for  diurnal  cycling  greatly  reduced  error  compared  to  using  a  constant  value  of  K^. 
Application  of  the  model  to  winter  open-water  and  ice-covered  periods  indicates  that  both  and 
are  higher  than  previously  thought.  Results  show  is  not  negligible,  and  probably  lies  between 

about  0.5  and  0.75,  while  is  on  the  order  of  0.1,  not  substantially  different  than  values  used  by 
others. 

The  results  of  this  study  indicate  that  photosynthetic  activity  can  be  important  in  winter  water 
quality.  The  presence  of  an  ice  cover  does  not  prevent  photosynthetic  activity.  These  data  support 
the  concept  that  improved  understanding  of  the  interactions  among  frazil  ice,  ice-covered  water, 
and  the  microbiology  of  ice  and  water  may  suggest  ways  to  improve  winter  water  quality.  The  use 
of  phospholipid  ester-linked  fatty  acid  analysis  of  samples  extracted  directly  from  water  may  avoid 
problems  with  culturing.  Water  samples  should  be  obtained  within  and  below  the  frazil  accumula¬ 
tion  in  order  to  more  clearly  identify  potential  community  differences.  Additional  field  data  are 
required  to  test  and  refine  the  model,  including  the  effects  of  air  temperature,  and  to  develop  a  range 
of  reaeration  and  deoxygenation  coefficients  useful  to  modelers. 

ACKNOWLEDGEMENTS:  The  authors  gratefully  acknowledge  Mr.  Rocky  Martin,  Public 
Works  Director,  Town  of  Hinesburg,  VT,  for  his  assistance  in  locating  and  installing  the  field 
sampling  site,  and  Ms.  Susan  E.  Hardy,  Environmental  Science  Branch,  CRREL,  for  her  input  and 
suggestions  to  this  technical  note. 

POINTS  OF  CONTACT:  This  technical  note  was  written  by  Ms.  Kathleen  D.  White,  Dr.  Charles 
M.  Reynolds,  Mr.  David  Ringelberg,  Ms.  Karen  L.  Foley,  and  Mr.  Lawrence  B.  Perry  of  the  U.S. 
Army  Engineer  Research  and  Development  Center’s  Cold  Regions  Research  and  Engineering 
Laboratory,  and  Dr.  Jeffrey  P.  Laible  of  the  University  of  Vermont.  For  further  information,  contact 
Ms.  White  (603-646-4187,  Kathleen.D.White@erdc.usace.army.mil)  or  the  managers  of  the  Water 
Operations  Technical  Support  Program,  Dr.  John  W.  Barko  (601-634-3654,  John.W.Barko® erdc. 
usace.army.mil),  or  Mr.  Robert  C.  Gunkel  (601-634-3722,  Robert.C.Gunkel@erdc.usace.army. 
mil).  This  technical  note  should  be  cited  as  follows: 

White,  K.  D.,  Reynolds,  C.  M.,  Laible,  J.  P.,  Foley,  K.  L.,  and  Perry,  L.  B.  (2001). 
“Low-temperature  microbial  activity  in  river  systems,”  Water  Quality  Technical  Notes 


19 


Water  Quality  Technical  Note  AM-1 1 
September  2001 

Collection  (WQTN  AM- 11),  U.S.  Army  Engineer  Research  and  Development  Center, 
Vicksburg,  MS.  www.wes.army.mil/el/elpubs/wqtncont.html 

REFERENCES: 

Chambers,  P.  A.,  Scrimgeour,  G.  J.,  Pietroniro,  A.,  Culp,  J.  M.,  and  Loughran,  I.  (1993).  “Oxygen  modelling  under  river 
ice  covers.”  Proceedings,  Workshop  on  the  Environmental  Aspects  of  River  Ice,  NHRI  Symposium  Series  No.  12, 
National  Hydrology  Research  Institute:  Saskatoon,  Saskatchewan,  T.  D.  Prowse,  ed.,  235-260. 

Cole,  T.  M.,  and  Buchak,  E.  M.  (1995).  “CE-QUAL-W2:  A  two-dimensional,  laterally  averaged,  hydrodynamic  and 
water  quality  model.  Version  2.0,”  Instruction  Report  EL-95-1,  U.S.  Army  Engineer  Waterways  Experiment 
Station,  Vicksburg,  MS. 

Ehcart,  J.  W.,  and  Park,  H.  (1989).  “Effects  of  temperature  variation  on  critical  stream  dissolved  oxygen,”  Water 
Resources  Research,  25(2),  1 45-151 . 

Fletcher,  M.  (1979).  “The  aquatic  environment.”  Microbial  Ecology:  A  Conceptual  Approach  (J.M.  Lynch  and  N.J. 
Poole,  eds.)  John  Wiley  and  Sons,  New  York,  92-1 14. 

Freeze,  R.  A.,  and  Cherry,  J.  A.  (1979).  Groundwater.  Prentice-Hall,  Inc.,  Englewood  Cliffs,  NJ. 

Gamble,  R.  W.  (1971).  Water  quality  on  the  St.  John  River  under  winter  conditions.  New  Brunswick  Department  of 
Fisheries  and  Environment,  Fredericton,  NB. 

Guasch,  H.,  Armengal,  J.,  Marti,  E.,  and  Sabater,  S.  (1998).  “Diurnal  variation  in  dissolved  oxygen  and  carbon  dioxide 
in  two  low-order  streams,”  Water  Research,  32  (4),  1067-1074. 

Hynes,  H.  B.  N.  (1970).  The  ecology  of  running  waters.  University  of  Toronto  Press,  Toronto. 

James,  W.  F.,  Eakin,  H.  L.,  Gunnison,  D.,  and  Barko,  J.  W.  (1992).  “Sediment-overlying  water  relationships  affecting 
wintertime  dissolved  oxygen  conditions  in  the  Big  Eau  Plaine  Reservoir,  Wisconsin.”  Water  Quality  Research 
Program  Technical  Report  W-92-2,  U.S.  Army  Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 

Landine,  R.  C.  ( 1 970).  “Prediction  of  dissolved  oxygen  levels  in  the  South  Saskatchewan  River  in  winter.”  Proceedings, 
International  Symposium  on  Water  Pollution  Control  in  Cold  Climates,  Water  Pollution  Control  Research  Series, 
U.S.  Environmental  Protection  Agency,  Washington,  DC,  96-112. 

Macdonald,  G.,  Holley,  E.  R.,  and  Goudey,  S.  (1989).  “Athabasca  River  winter  reaeration  investigation.”  Report 
prepared  for  the  Environmental  Assessment  Division,  Alberta  Environment,  Edmonton,  AB. 

Marchand,  P.  J.  (1996).  “Life  in  the  cold:  An  introduction  to  winter  ecology,”  University  Press  of  New  England,  Hanover, 
NH. 

McBean,  E.,  Farquhar,  G.,  Kouwen,  N.,  and  Dubek,  O.  (1979).  “Prediction  of  ice-cover  development  in  streams  and  its 
effect  on  dissolved  oxygen  modeling,”  Canadian  Journal  of  Civil  Engineering  6,  197-207. 

Pietroniro,  A.,  Chambers,  P.  A.,  and  Ferguson,  M.  E.  ( 1 998).  “Application  of  a  dissolved  oxygen  model  to  an  ice-covered 
river,”  Canadian  Water  Resources  Journal  23  (4)  351-368. 

Reitncr,  B.,  Herzig,  A.,  and  Herndl,  G.  J.  (1997).  “Microbial  activity  under  the  ice  cover  of  the  shallow  Neusiedler  See 
(Austria,  Central  Europe),”  Hydrobiologia  357,  173-184. 

Reynolds,  C.  M.,  Perry,  L.  B.,  Pidgeon,  C.  S.,  Martel,  C.  J.,  and  Koenen,  B.  A.  (1998).  “Effects  of  diurnal  temperature 
cycles  on  carbon  dioxide  evolution  from  soil.”  49th  Arctic  Science  Conference  [25-28  October,  1998,  Fairbanks, 
AK]. 

SAS  Institute,  Inc.  (1996).  “SAS®  SoftwareVersion  6.12,”  Cary,  NC. 

Sasser,  M.,  and  Wichman,  M.  D.  (1991).  “Identification  of  microorganisms  through  use  of  gas  chromatogoraphy  and 
high-performance  liquid  chromatography.”  Manual  of  Clinical  Microbiology,  5th  ed.  W.J.  Hausler,  Jr.,  K.L. 
Herrmann,  H.D.  Isenberg,  H.J.  Shadomy,  ed.,  American  Society  for  Microbiology,  Washington,  DC,  111-118. 

Schreicr,  H.,  Erlebach,  W.,  and  Albright,  L.  (1980).  “Variations  in  water  quality  during  winter  in  two  Yukon  rivers  with 
emphasis  on  dissolved  oxygen  concentration,”  Water  Research,  14,  1345-1351. 


20 


Water  Quality  Technical  Note  AM-1 1 
September  2001 


Shallock,  E.  W.,  and  Lotspeich,  F.  B.  (1974).  Low  winter  dissolved  oxygen  in  some  Alaskan  rivers.  U.S.  Environmental 
Protection  Agency  Arctic  Environmental  Research  Laboratory:  College,  AK. 

Simonsen,  J.  F.,  and  Harremoes,  P.  (1978).  “Oxygen  and  pH  fluctuation  in  rivers.”  Water  Research,  12, 477-489. 

Streeter,  H.  W.,  and  Phelps,  E.  B.  (1925).  “A  study  of  the  pollution  and  natural  purification  of  the  Ohio  River,  III.  Factors 
concerned  in  the  phenomena  of  oxidation  and  reaeration,”  United  States  Public  Health  Service  Public  Health 
Bulletin  No.  146,  U.S.  Government  Printing  Office:  Washington,  DC. 

Tchobanoglous,  G.,  and  Schroeder,  E.  D.  (1987).  Water  quality.  Addison-Wesley  Publishing  Co.,  Inc.,  Reading,  MA. 

Thomann,  R.  V.,  and  Mueller,  J.  A.  (1987).  Principles  of  surface  water  quality  modeling  and  control.  Harper  Collins 
Publishers,  Inc.,  New  York. 

U.S.  Environmental  Protection  Agency.  (1995).  “QUAL2E— The  Enhanced  Stream  Water  Quality  Model.”  U.S. 
Environmental  Protection  Agency  Software  Package  EPA-823-C-95-006,  Documentation  EPA-823-B-95-003, 
Center  for  Exposure  Assessment  Modeling,  Athens,  GA. 

Vincent,  W.  F.,  Rae,  R.,  Laurion,  I.,  Howard-Williams,  C.,  and  Priscu,  J.  C.  (1998).  “Transparency  of  Antarctic 
ice-covered  lakes  to  solar  UV  radiation,”  Limnology  and  Oceanography,  43(4),  618-624. 

White,  K.  D.  (1998).  “1-D  streamwise  finite  element  model  of  dissolved  oxygen  under  river  ice.”  Ice  in  Surface  Waters, 
Proceedings,  14th  International  Symposium  on  Ice.  July  27-31,  1998,  Potsdam,  NY. 

White,  K.  D.,  and  Melloh,  R.  A.  (1999).  “Diurnal  variation  in  dissolved  oxygen  measurements  during  late  winter 
ice-covered  period.  Sleeper’s  River,  VT.”  Proceedings,  10th  International  Conference  on  Cold  Regions  Engineer¬ 
ing,  August  16-19,  Lincoln,  NH,  550-560. 

Whitfield,  P.  H.,  and  McNaughton,  B.  (1986).  “Dissolved-oxygen  depressions  under  ice  cover  in  two  Yukon  rivers,” 
Water  Resources  Research  22  (12),  1675-1679. 


NOTE:  The  contents  of  this  technical  note  are  not  to  be  used  for  advertising,  publication, 
or  promotional  purposes.  Citation  of  trade  names  does  not  constitute  an  official  endorse¬ 
ment  or  approval  of  the  use  of  such  products. 


21 


20010928  006 


96  qnr 

31310590  30V  SNOtllOg  SnOiA3Hd  _  lujoj  OtlQ 


z 

lU 

•as^aX^-^  ^^5  paAOJcIdy 

iN3W3iViS  NOlinaiHiSIQ  ’Z 

( 

uaawnN  j.ovhinoo  agoNn  aaavdaud  a 

V  s 

'} 

ifl  > 

HHawnN  laodaa  aoiiNOW  o 

JU 

IT-W\"  330N  IBOT^qosi 

WV^  X 

UaBHON  HO/QNV  3nili  iHOdSa  -a 

00  7. 

iBjuBO  ABsi  >s  tfoaBBse^  aaBUxgua  Auay  *S*I1 

uo 

AON3DV  DNiiVNlOiaO  V 

NOlivwyodNi  oNiAd'iNaai  iBOday  ‘v 

eie 

d  ^  I J-'-'  ■* 


OliO 


T0-<7Z-6 

3iva 


u3awnN  av 


